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A MONTHLY VECTOR WIND PROFILE MODEL 


1. INTRODUCTION 

The purpose of this user’s guide is to provide an understanding of the background, theoretical 
concepts and methodology for construction of vector wind profiles based on a statistical model. 
With this understanding the user can embark on the intended application, to provide the launch 
vehicle engineering design community a product that has wide application in establishing realistic 
estimates of the dispersions of critical vehicle design parameters related to wind profile 
dispersions. This user’s guide includes a description of the computational process for the model 
including specification of input data, parameter settings, and output data formats. Sample output 
data listings are provided to aid the user in the verification of test output generated by the user. 


2. BACKGROUND 

The most useful engineering design application of a wind profile model is the establishment of 
preliminary design ranges for angle of attack, a, angle of sideslip, p, aerodynamic pressure, q, and 
the two aerodynamic load indicators, products, qoc, and qp. These and other flight variables are 
derived from ascent flight 6-degree of freedom trajectory simulations using wind model profiles. 
The trajectory variables are used in the evaluation of load indicators at locations of expected 
vehicle wind sensitivity. A load indicator is an algorithm that relates an external loads such as q to 
stress at a specific point on the vehicle structure; for the Space Shuttle, the algorithms are for 
rigid body loads. 1 Elastic body loads are determined from flutter and vibration analyses using 
model wind profiles augmented to include small scale wind perturbations. Another useful 
application is the estimation of flight performance reserve (FPR) for propellant to ensure orbital 
insertion by protecting for flight dispersions attributable in part to wind profile dispersions. ’ 1 
Following the preliminary vehicle design using a wind profile model, trade studies arc made to 
establish a requirement to bias steering to reduce wind loads. The usual procedure is to establish 
first-stage steering based on the profile of monthly mean winds in the pitch and yaw planes. Wind 
profile models were developed for alternatives other than the monthly means. 3 When sufficient 
engineering data are established, structural loads and performance assessments are made using 
samples of high resolution wind profile measurements. Currently, for Kennedy Space Center 
(KSC) this data sample is 150 Jimsphere profiles per month. 


3. APPLICATION 

The monthly vector wind profile (VWP) model is suitable for applications in preliminary launch 
vehicle design studies that require assessments of vehicle trajectory and aerodynamic loads 
dispersions attributable to monthly wind profile dispersions. Launch vehicle ascent guidance and 
control system (auto-pilot) steering commands are programmed for flight through the profile of 
monthly mean wind. Various vehicle programs have used different terms for the programmed 
steering commands. The early Saturn program called the commands the chi-tilt program, later the 
wind-biased trajectory, and the Space Shuttle program calls the steering commands the I-Load. 



During initial studies a number of months are used to establish the model profiles that produce the 
largest monthly dispersions of ascent vehicle aerodynamic load indicators, Qa and Qp. Because 
the largest monthly dispersions for wind occur during the winter high-wind months, it is 
appropriate to use the worst month from the winter season for establishment of the design 
reference, Qa and Qp, dispersions. 


4. THEORETICAL CONCEPTS 


The theoretical basis for the vector wind model is the wind components at any two altitudes 
within the altitude range of the model are quadravariate normal. 

A quadravariate normal distribution for wind vectors, {U } , V]}, at an altitude, Zj, and wind 
vectors, {U 2 , V 2 }, at an altitude, Z 2 , has the property that the conditional wind vectors at an 
altitude, Z 2 , are bivariate normally distributed given specific values for the wind vector, 

{ u i- Vf}, at altitude, Z { . 


Symbolically, this is 


f (u 2 ,v 2 |u;,v;) 


r (U), Vj, u 2 , v 2 ) 
t- (u;. v;) 


( 1 ) 


The 14 statistical parameters for the quadravariate normal distribution for the zonal, U, and 
meridional, V, wind components are the four mean values, U,, Vi, U 2 , andV 2 ; the four standard 
deviations, SDUi, SDVi, SDU 2 , and SDV 2 ; and the six correlation coefficients for like and unlike 
variables between altitudes Zi and Z 2 : R(Ui,U 2 ), R(Vj,V 2 ) and R(Ui,V 2 ) and R(Vi,U 2 ), 
respectively, and unlike variables at the same altitude Zi or Z 2 :R(Ui,Vi) and R(U 2 ,V 2 ) . 


5. WIND PROFILE CONSTRUCTION 

The procedure is (1) define 12 specific wind vectors to the 99 percent probability ellipse at a fixed 
reference altitude, (2) compute the five parameters for the conditional bivariate normal probability 
distributions for all altitude levels above and below the reference altitude, and (3) find the 
intercept to the conditional probability ellipse toward the mean values at the reference altitude. 
This conditional wind vector closely approximates the largest vector wind shear between the 
reference altitude and each ol the other altitudes. 

For engineering applications, the conditional wind vectors arc expressed in polar coordinates as 
wind speed and wind direction in the standard meteorological coordinate system. Thus, the 
vector wind profile model is defined by the 12 equally spaced (30-degree) increments Irom the 
centroid of the 99 percent probability ellipses at each reference altitude, which are the given 
values for the 99 percent conditional wind vectors that yield the largest shear at all other altitudes 
above and below the reference altitude. 
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For the N altitudes for the available data base there are 12 by N vector wind profiles from the 
surface (station elevation) to the maximum altitude. These vector wind profiles as a function of 
altitude above mean sea level are expressed as wind speed and wind direction in the standard 
meteorological coordinate system. The vector wind profiles are derived from the vector wind 
profile model for a selected month, called the design reference month, representative of the high 
wind months. 

This coordinate system is chosen because a vehicle could have a flight azimuth different from 
either axis associated with the wind components. Furthermore, for operations the wind data is 
provided in the meteorological coordinate system. 

The wind coordinate system used in the Space Shuttle program as a function of flight azimuth, 
FA, is 

W X (FA) = W cos (0-FA) [ in-plane wind component] (2) 

where, a headwind is a positive in plane wind component and a tailwind is negative. 

W y (FA) = W sin (0 - FA) [out-of-plane wind component] (3) 

where, a right-to-left out-of-plane wind component is a positive crosswind and left-lo-right is 
negative crosswind. The wind direction, 0, measured in degrees clockwise from true north, is the 
direction from which the wind is blowing. W is wind speed and FA is flight azimuth measured in 
degrees clockwise from true north. 


6. EQUATIONS 

This section presents the specific equations to compute (1) the given wind vectors at the reference 
altitude, Z], and (2) the five conditional bivariate normal parameters at Z 2 , i.e., the conditional 
component means, the conditional component standard deviations, and the conditional correlation 
coefficients. 1 The five conditional statistical parameters are used to compute the conditional 
oivariate normal 99 percent probability ellipse from which the conditional wind vector that 
approximates the largest shear between the reference altitude, Z|, and all other altitudes Z 2 above 
and below the reference altitude. Notation used for the statistical parameters is the following: 
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Means 

Variances 

Standard 

Correlation 

Covariances 



Deviations 

Coefficients 

(m 2 /s 2 ) 

(m/s) 

(m 2 /s 2 ) 

(m/s) 

(unitless) 

MU, 

SU, 

SDU, 

RU,V, 

su,v, 

MV, 

sv, 

SDV, 

RU,V 2 

su,v 2 

mu 2 

su 2 

sdu 2 

rv,u 2 

sv,u 2 

mv 2 

SV 2 

sdv 2 

ru 2 v 2 

SU 2 V 2 




ru,u 2 

SU,U 2 




rv,v 2 

sv,v 2 


The general expression for the covariance is SUjVj = RU i Vj(SDU i )(SDVj) . 

The first step is to define the 12 given wind vectors to the 99 percent probability ellipse at the 
reference altitude. These given wind vectors from the centroid of the probability ellipse are 
defined at 30-degree increments of clocking angle, CA, measured in the standard mathematical 
convention (counterclockwise). 


where, 


u; = MU, + RS cos(CA) 

(4) 

v; = MV, + RS sin(CA) 

(5) 

It 

c sy 

& 


A V- 2 In (1 - P) 



where P is probability = 0.99 , and 


A 2 = 


1 

l-(RU,Vi) 2 


" cos(CA) f 

, SDUj J 


2RU,V, cos (CA) sin( CA) 
SDU, SDV, 


sin(CA) Y 

SDV, J 


( 6 ) 


The conditional mean vectors, {CMU 2 , CMV 2 }, at altitude, Z 2 , given specific wind 
vectors, { UpV,* }, at the reference altitude are 


CMU 2 =MU 2 + (T, + T 2 ) / [1 - RU,V, * RU 1 V 1 )] 


(7) 


where 

T, = [RU,U 2 - RU, V 2 * RU,Vi] * ( uj - MU,) (SDU 2 / SDU,) 


and 

t 2 = [ru,v 2 - ru,u 2 * ru,v,i * (v; - mv,) (sdu 2 / sdv,) 

CMV 2 = MV 2 + (T 3 + T 4 ) / [ 1 - RU,V, * RU,V,] («) 

where, 

T 3 = [RV,U 2 - RV,V 2 * RU 1 V 1 ] * (u;-MU.) * (SDV 2 /SDU,) 

and 
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T 4 = [RV,V 2 - RV 1 U 2 * RU 1 V 1 ] * (V,* - MV,) * (SDV 2 / SDV,) 


The conditional standard deviations are 


CSDU 2 = [sigma (1,1)] 1/2 

(9) 

and 


CSDV 2 = [sigma (2,2)] 1/2 

(10) 


The conditional correlation coefficient is 

CRU 2 V 2 = sigma (1,2)/ (CSDU 2 )(CSDV 2 ) (11) 


where 

sigma (1,1) = SU 2 - SU,U 2 [SU,U 2 * SV, - SU 2 V, * SU,V,] / D 

- SV,U 2 [ - SU,U 2 * SU, V, + SU 2 V, * SU,] / D 

sigma (2,2) = SV 2 - SU,V 2 [SU, V 2 * SV, - SV,V 2 * SU,V,] / D 

- SV,V 2 [ - SU,V 2 * SV,U 2 + SV,V 2 * SU,] / D 

sigma (1,2) = SU 2 V 2 - SU,V 2 [SU,U 2 * SV, - SU 2 V, * SU,V,] / D 

- SV,V 2 [SU,U 2 * SU,V, + SU 2 V, * SU,](1/D) 

and, D = (SU,)(SV,) - (SU|V,)(SU,V,). 

Note that the given wind values U[ and V 2 are required for the conditional mean component and 
not for the conditional standard deviations. 

A vector wind model profile consists of the given wind vector at Z, and wind vectors UC, and 
UC 2 calculated from each conditional 99 percent probability ellipse (PE) at all other altitudes Z 2 
above and below Zl. Selection of UC1 and UC 2 on the conditional PE at a clocking angle 180 
degrees from the clocking angle of the given wind vector produces the desired near maximum 
vector shear between Z, and Z 2 . 

UC 2 = CMU 2 + RSC cos(CC) ( 1 2) 

and 

VC 2 = CMV 2 + RSC sin(CC) ( 1 3) 


where 

CC = CA + 180, i.e., the clocking angle to the given vector plus 180 degrees 

RSC = p=J= = 

A c V-2 1n (1-P) 
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where 


P is probability = 0.99, and, 


A c 2 = 


cos(CC) ) 2 2CRU 2 V 2 cos(CC) sin(CC) f sin(CC) f 


1-(CRU 2 V 2 ) 2 


^ csdu 2 j 


csdu 2 *csdv 2 


l csdv 2 J 


(14) 


The schematic in figurel illustrates the construction of a given wind vector at 12-km altitude for a 
clocking angle of 30 degrees and the conditional wind vector at 10-km, which is at a clocking 
angle of 210 degrees, (30 + 1 80), on the 99 percent conditional probability ellipses. 



I50.0 -25.0 0.0 25.0 SO.O 75,0 100.0 125.0 

U, m/s 


Figure 1. Schematic of Profile Construction Between a Reference Altitude of 12 km and an Altitude 
of 10 km, Clocking Angle 30°, 

The monthly mean profile and the 12 wind profiles for KSC, February for a reference altitude of 
12-km are illustrated in figure 2 for the zonal wind component (U) and figure 3 for the meridional 
wind component (V). 



altitude, km 



Figure 2. Mean Profile (thick solid) and Vector Wind Model Profiles for Z 0 = 12 km, 

U-Component KSC, February 



Figure 3. Mean Profile (thick solid) and Vector Wind Model Profiles for Z 0 = 12 km, 

V-Component KSC, February 
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7. COMPUTATIONAL PROCESS 


Preferably, the input data consist of a statistically representative sample of twice daily wind 
profiles obtained during at least a ten-year period of record yielding approximately 600 profiles 
per month. The database used for KSC in the example presented in this guide is the 19-year 
(1956-74) serially complete Rawinsonde profile set established by NASA for the Space Shuttle 
Program. 

The computer code for the Monthly Vector Wind Profile Model consists of three subroutines. 
The first subroutine, QUAD, uses the input data file for calculations of the monthly quadravariate 
normal statistics at any reference altitude, HO, and any other altitude, H, above or below HO, 
within the altitude range of the available data base. The quadravariate normal statistics are used in 
the second subroutine, CSTAT, for the calculations of the conditional bivariate normal statistics at 
H given a wind vector to the pth percentile ellipse at HO. The 12 given wind vectors at each HO 
and the conditional bivariate normal statistics are used in the third subroutine, PRO, lor the 
construction of 12 vector wind profiles for each HO. 


Requirements and specifications for input parameters and input and output data files are included 
in the process description below. The computer codes for the various subroutines are not 
included in this report. 

The driver for the computations contains dimension statements for the derived statistical arrays of 
the form illustrated below for 28 altitude levels. 

DIMENSION XB (28) ,YB(28) ,SDX{28) ,SDY(28) ,RXY(28) ,DET1(28,28) ,DET2(28,28) 

C 

C SET MONTH AND NUMBER OF LEVELS 
C IMM=MONTH , 

C NLV=NUMBER OF ALTITUDE LEVELS FROM SURFACE TO MAXIMUM ALTITUDE USUALLY BUT NOT RESTRICTED 
C TO 1KM INTERVALS 

C 

DATA IMM/2/ 

DATA NLV/28/ 

C 

C CALCULATE QUADRAVARIATE NORMAL STATISTICS 

CALL QUAD (XB, SDX, YB, SDY, RXY , DET1 , DET2 , NLV, IMM) 

C 

C CALCULATE CONDITIONAL STATISTICS 

CALL C ST AT ( XB , SDX , YB , SDY , RXY , DET 1 , DET2 , NLV, IMM ) 

C 

C CONSTRUCT PROFILES 

CALL PRO {NLV, IMM) 

STOP 

END 

C 

SUBROUT I NE QUAD ( XB , SDX , YB , SDY , RXY , DET 1 , DET 2 , NLV , I MM ) 

C THIS PROGRAM GENERATES THE QUADRAVARIATE NORMAL STATISTICS REQUIRED 
C FOR COMPLETE VECTOR WIND PROFILE MODEL FOR A SELECTED MONTH 

C 

C NOTE: DIMENSION STATEMENTS BELOW ARE FOR A MAXIMUM OF 28 LEVELS , WHICH IS THE 
C PRACTICAL LIMIT FOR RAWINSONDE PROFILES. 

C 

DIMENSION SX { 2 8 ) , SY{28) , SSX(28) , SSY (28) , SSXY (28) , XB(28) , YB(28) , 

$SDX { 28 ) , SDY (28) ,RXY{28) ,SXXP{28,28) ,SYYP{28,28) ,SXYP(28,28) , 

$SYXP(28,28) , RXXP<28, 28) ,RYYP{28,28) ,RXYP(28, 28) ,RYXP{28, 28) , 

$ IDD < 5 ) , IWD (28), IWS (28) , X ( 2 8 ) , Y { 2 8 ) , KK ( 2 8 ) , DETl (28,28) , DET 2 (28,28) 

CHARACTER* 1 0 MON , MONTH ( 1 2 ) 
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DATA MONTH/ 'JANUARY ' , ‘FEBRUARY 1 , 'MARCH ' , 'APRIL 

$ ‘ MAY ' , ' JUNE * , ' JULY ' , ' AUGUST ' , * SEPTEMBER 1 , 

$ ' OCTOBER * , ' NOVEMBER ' , 1 DECEMBER • / 

C 

C PARAMETERS TO BE SPECIFIED IN DATA STATEMENTS BELOW 
C SFC=ELEVATION IN KM TO 3 DECIMAL PLACES (IE. TO NEAREST METER) 

C IOUTl=l FOR BIVARIATE NORMAL STATISTICS OUTPUT TO UNIT 97 
C IOUT2=l FOR INTERLEVEL CORRELATION COEFFICIENTS BETWEEN LIKE COMPONENTS 
C OUTPUT TO UNIT 98 

C IOUT3 = 1 FOR INTERLEVEL CORRELATION COEFFICIENTS BETWEEN UNLIKE COMPONENT 
C OUTPUT TO UNIT 99 

C 

DATA SFC/Q . 003 / 

DATA IOUT1/1/ 

DATA IOUT2/1/ 

DATA IOUT3/1/ 

C 

C INPUT DATA FROM UNIT 2 IS SERIALLY COMPLETE RAWINSONDE DATA BASE 
C FORMAT OF INPUT DATA IS GIVEN IN STATEMENT 450 BELOW 

C INPUT PILE 

OPEN (UNIT»2 , FILE* ' SERCOMKSC . TXT ' ) 

C OUTPUT PILES (IP NEEDED) 

C NOTE •2" IN OUTPUT FILE NAME SPECIFIES THE APPLICABLE MONTH , FEBRUARY FOR THIS EXAMPLE 
OPEN (UNIT* 9 7, FI LE= 'BNSKSC2.TXT ' ) 

OPEN ( UNIT* 9 8 # PI LE» • CORAKSC2 . TXT 1 ) 

OPEN { UNIT= 9 9 , FILE* 1 CORBKSC2 . TXT • ) 

C 

SUBROUTINE CSTAT (XB, SDX , YB , SDY, RXY, DET1, DET2 , NLV, IMM) 

C 

C CALCULATES CONDITIONAL BIVARIATE NORMAL STATISTICS FOR PTH PERCENTILE 
C CONDITIONAL ELLIPSES AT ALL ALTITUDES H r FOR EACH OP 12 GIVEN VECTORS 
C TO THE PTH PERCENTILE ELLIPSE AT REFERENCE ALTITUDE HO (ALL ALTITUDES 
C SURPACE TO MAXIMUM ALTITUDE) . NOTE THAT P IS SET BELOW. 

C 

DIMENSION XB (28) ( SDX(28) ,RXY(28) , YB(28) ,SDY(28) ,DET1{28,28) , 

*DET2(28 / 28} ,RYYPUP(28) ,RXXPDN(28) ,RXXPUP(28) ,RYYPDN(28) , RYXPUP(28) 

* , RXYPDN (28) , RXYPUP (28 ) ,RYXPDN{28) ,XGG(28) ,YGG{28) , YYCLK (28, 28) , 

*CBNUP ( 28 , 28,5) ,CBNDN(28, 28, 5) ,XCLK(28) ,YCLK(28) , XXCLK ( 28 , 28 ) , 

*RCC (12) 

CHARACTER IMO*10 , ID1 *7 , ID6*15 
CHARACTER* 10 IMONN(12) 

DATA lUOm/ 'JANUARY ‘ , 'FEBRUARY 1 , 'MARCH ' , 'APRIL 
* ' MAY 1 , 1 JUNE ’ , ' JULY T , ' AUGUST ' , ' SEPTEMBER ' , 

* ’ OCTOBER ' , ' NOVEMBER ' , ’ DECEMBER ' / 

DATA ID1 / ' KSC ' / 

DATA ID6/' GENERAL CASE '/ 

C SET PROBABILITY LEVEL FOR DERIVATION OF GIVEN WIND VECTOR TO PTH PERCENTILE 
C ELLIPSE AT HO 

DATA P/.99/ 

DATA IOUT4/1/ 

C OUTPUT FILE IF NEEDED SET IOUT4=l 

C NOTE "2" IN OUTPUT FILE NAME SPECIFIES THE APPLICABLE MONTH 
OPEN (UNIT=27 , FILE* 1 CSTATKSC2 . TXT 1 ) 


SUBROUTINE CONDELL (XBO , YBO , SDXO , SDYO , RXYO , XB2 , YB2 , SDX2 , SDY2 , RXY2 , 

* RXXP , RYYP , RXYP , RYXP , XG , YG , CMX , CMY , C SDX , CSDY , CRXY ) 

DIMENSION SIGMA (2,2) , XU ( 2 ) ,VX(2) , US < 2 ) ,COR(2, 2) , VI (2,2) , 

*V2 (2, 2) , V3 (2, 2) , V4 (2, 2} , Ul (2) , U2 (2) ,X{2) , W24 (2,2) ,VW{2, 2) ,SD(4) , 
*RHO (4,4) , U { 4 ) , S (4 , 4 ) 






C DEFINITION OF INPUT VARIABLES 

C U ( 1 ) =MEAN OF X COMPONENT AT ALT HO 
U(1)=XB0 

C U ( 2 ) =MEAN OF Y COMPONENT AT ALT HO 
U ( 2 ) = YB0 

C U ( 3 ) =MEAN OF X COMPONENT AT ALT H2 
U ( 3 ) =XB2 

C U (4 ) =MEAN OF Y COMPONENT AT ALT H2 
U (4 ) =YB2 

C SD { 1 ) =STD . DEV . OF X AT ALT HO 
SD { 1 ) =SDX0 

C SD ( 2 ) =STD . DEV . OF Y AT ALT HO 
SD (2) =SDY0 

C SD { 3 ) =STD . DEV . OF X AT ALT H2 
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SD ( 3 ) =SDX2 

C SD ( 4 ) =STD . DEV . OF Y AT ALT H2 
SD{4)=SDY2 

C RHO (1,2) CORRELATION OF X AT HO WITH Y AT HO 
RHO < 1,2) =RXY0 

C RHO ( 1 , 3 ) CORRELATION OF X AT HO WITH X AT H2 
RHO { 1 , 3 ) =RXXP 

C RHO ( 1 , 4 ) CORRELATION OF X AT HO WITH Y AT H2 
RHO (1,4) =RXYP 

C RHO ( 2 f 3 ) CORRELATION OF Y AT HO WITH X AT H2 
RHO {2,3} =RYXP 

C RHO (2 f 4} CORRELATION OF Y AT HO WITH Y AT H2 
RHO { 2 , 4 ) = RYYP 

C RHO (3,4) CORRELATION OF X AT H2 WITH WITH Y AT H2 
RHO (3,4) =RXY2 

C X { 1 } =X COMPONENT OF THE GIVEN VECTOR 
X { 1 ) =XG 

C X ( 2 } =Y COMPONENT OF THE GIVEN VECTOR 
X { 2 } =YG 

C DEFINITION OF OUTPUT VARIABLES 
C US (1) CONDITIONAL MEAN OF X COMPONENT 
C US { 2 ) CONDITIONAL MEAN OF Y COMPONENT 

C COR (1,1) CONDITIONAL STANDARD DEVIATION OF X COMPONENT 
C COR (2,2) CONDITIONAL STANDARD DEVIATION OF Y COMPONENT 
C COR (1,2) CONDITIONAL CORRELATION OF X AND Y COMPONENT 

Q* ********** ************** ******************************** 

SUBROUTINE PRO{NLV,IMM) 

C MONTHLY WIND PROFILE MODEL 

C CLOCKING ANGLE OF CONDITIONAL WIND VECTOR IS 180 DEG FROM CLOCKING 
C ANGLE OF GIVEN WIND VECTOR 

DIMENSION XC{28) ,YC{28) , VWS (28),XB(28),YB(28), SDX (28), 

*5DY (2 8) , RXY (28) , SPWS<28) ,SPWD(28) , ICLA(28) , 

*SPWSTE{336, 28) , SPWDTE { 3 3 6 , 28 ) , 

*RPYD( 115000) ,RPYU (115000) , 

* RPXD (1150 0 0) , RPXU (115000) , DATWD{28) , DATWS(28) , 

*DATWDD {33 6 , 28 ) , DATWSS (336,28) , UU (33 6 , 28 } ,W(336,28) 

C SET PROBABILITY LEVEL FOR CONDITIONAL ELLIPSE 
DATA P/.99/ 

CHARACTER IMO*10 
CHARACTER* 10 IMONN(12) , IMON 

DATA IMONN/ 'JANUARY 1 , 'FEBRUARY ’ , 'MARCH ' , 'APRIL 

* - MAY 1 , ' JUNE ' , ' JULY ' , ' AUGUST ‘ , * SEPTEMBER 1 , 

* ' OCTOBER ' ! ' NOVEMBER ' , ' DECEMBER ’ / 

CHARACTER IDl*7,ID5*15 
CHARACTER IDFI*3 , IH0I*3 
DATA IDl/'KSC '/ 

DATA IDS/ ' '/ 

C INPUT GIVEN WIND VECTORS AND CONDITIONAL BIVARIATE NORMAL STATISTICS 
OPEN (UNIT=27 , FILE= * CSTATKSC2 . TXT ' ) 

C OUTPUT FILE: UNIT 86 LISTS WIND PROFILES AND SHEARS FOR EACH REFERENCE ALTITUDE HO 
C (0 TO (NLV-l)KM) WITH RESPECT TO 12 CLOCKING ANGLES (0 TO 330 DEG AT 30 DEG INCREMENTS) 
C OPEN {UNIT** 8 6 , FILE* 1 PROKSC2 . TXT * ) 

C OUTPUT FILE UNIT 72 LISTS WIND PROFILES IN SPREAD SHEET 673 COLUMN FORMAT 


C 
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CONTENTS 










C 

1 
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c 

50-73 
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= 2 KM 










C ETC. 

C (HO+1) *24-22 TO (HO+1) *24+1 SAME AS 2-25 EXCEPT H0=H0 
C 650-673 SAME AS 2-25 EXCEPT HO = 27KM 


II 



8. SAMPLE OUTPUT DATA 


Samples of output data from the three subroutines are presented in tables 1 through 5. 


Table 1 Bivariate Normal Statistics 

Source: Subroutine QUAD 

Input File: SERCOMKSC.TXT (KSC serially complete Rawinsonde, 1956-74) 
Output File: BNSKSC2.TXT 

BIVARIATE NORMAL STATISTICS FEBRUARY 
KSC SC RAWINSONDE, 1956-74 
NUMBER OF PROFILES: 1074 


z 

MEAN U 

S.D. U 

R (U, V) 

MEAN V 

S.D. V 

KM 

M/S 

M/S 


M/S 

M/S 

003 

0 . 653 

3.274 

-0.261 

-0.213 

3 . 676 

1 

3 . 699 

7.167 

-0.028 

1.620 

6.793 

2 

7.880 

7.741 

0 . 008 

1.484 

6 . 821 

3 

11.704 

8.195 

0 . 044 

1.684 

7.400 

4 

15.212 

9 . 109 

0 . 034 

2 . 105 

8 . 095 

5 

18 . 972 

10.157 

0.050 

2 .487 

8.979 

6 

22 . 948 

11.130 

0.112 

2.979 

9.612 

7 

26 . 584 

12.427 

0.165 

3.307 

10.478 

8 

30 . 235 

13 .789 

0 . 199 

3 .398 

11.539 

9 

34.256 

15.369 

0.215 

3 .500 

12 . 630 

10 

38 . 184 

16.446 

0.216 

3 .387 

13 . 883 
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42.127 

17.071 

0.228 

3 .324 

14 .900 

12 

44 . 840 

16 . 526 

0.227 

3 .486 

14 . 555 

13 

44.763 

15 . 053 

0.286 

3 . 525 

12 . 840 

14 

41 . 654 

13 . 073 

0.276 

3.330 

11 . 048 

15 

36,734 

11.440 

0.206 

2.902 

9.279 

16 

31.588 

10.264 

0.149 

2.600 

8.243 

17 

25.364 

9 . 197 

0.143 

1 . 939 

7.038 

18 

18.777 

8.485 

0.238 

1.407 

5 . 667 

19 

12.767 

7.840 

0.228 

0.989 

4.522 

20 

7 . 849 

7 .395 

0.254 

0.629 

3.887 

21 

5.211 

7 .254 

0.232 

0.179 

4.230 

22 

4 . 039 

7 . 653 

0.234 

-0 . 140 

4.103 

23 

3 .465 

7 . 862 

0.274 

-0.023 

4 . 098 

24 

3 .646 

8.261 

0.280 

0 . 091 

3.891 

25 

3 .881 

9 . 142 

0.347 

-0 . 020 

3.846 

26 

4.478 

9.819 

0.307 

0.111 

4 . 091 

27 

5 . 140 

10.570 

0.230 

0.352 

4.132 
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Table 2 Correlation Coefficients Between Like Wind Components (Page 1 of 2) 
Source: Subroutine QUAD 

Input File: SERCOMKSC.TXT (KSC serially complete Ra win sonde, 1956-74) 
Output File: CORAKSC2.TXT 

INTERLEVEL CORRELATION COEFFICIENTS BETWEEN LIKE WIND COMPONENTS 


Ho^nM/iHooTfcon^O'jHocooooffirtwr'MyjM^ioco 

oonrt^uim^^r'hoocooiocOhhLnmTiinniNfNoq^H 


OOOOOOOOOOOOOOr— IOOOOOOOOOOOOO 

MOrlTj^n(Ti(NMn(NHl/lr N DhrOU)rCO(NHOOinM^rl(J\ l Xl 

Hmr , iHi/iMcrir'»jmuioir)HOhrifOini/io\ocOhH(Ncoi^^ 

OO(NM^^inU)U)r'C 0 C 0 ^OC 0 C 0 hlJ>l/l^^(N(N(N(NHHH 


oooooooo O O O O O i — i O o o o o o o o o o o o o o 


fNfNHfO'^(NrOr-rOU)(^MCOO[^^l^<?\n^U)(N<ytVD^^OVDr' 

tHor^-cocxcctnrHootNt^r-roo^o^rof^ooTrinmnvo^^Ofof^ 

OOCNr0’<tfinVDVDr-r'000\OC0r'rVDV0in*3‘C0CNCNi-Hr-lrHrHrH 


OOOOOOOOOOOOtHOOOOOOOOOOOOOOO 


rH<j\a\r-fNt^inoor'<yvTj*uDOLnooo^rHr^cocof^(No>o^o\<y\ 

HOhVOHO\h^rl^U)^OrO'JLnCO(NUiy 5 rl(NHOiriO\^)ron 

OOCN^Tjtnvo^[sa 3 o^O(TlCOr'W^Ln^^^(NnHHHHH 


OOOOOOOOOOOHOOOOOOOOOOOOOOOO 


o(Ni/iLn^Hrnh^cor-oi/)oohfor'inhrivciNh^noooai 

T-imtHTHUDLnrnor''^roo'^coorH^Doro^J | OHchr'(Nr-Lnroro 


OOOOOOOOOOrHOOOOOOOOOOOOOOOOO 


<Tcor-- , 'ProinHCNrHoocNvoor'Coo>cNcoovcoinr''tnrHinr'Om 

COVDr-rOCNrHO\VO^O'<tf[-'CNincOCNCOCNcO<TOr^TpO\'<tfCNCNH 

oHfoinvor-r'Oo^o^cocoryjvoiMrt^MMHHOHHtHH 


OOOOOOOOOrHOOOOOOOOOOOOOOOOOO 


CO ’'P 

<r 

CN 

00 

CN 

rH 

CN 

10 

o 

o 

^p in 

rH 

^P 

r- 

o r- r- 

rH 


CO 

in 


[■" 

CN 

O 

VD 

O 


CN 

^P 

o\ 

<T 

co 

10 

CO 

o 

in 

r- 

O 

ID 

rH 

"P 

o in o 

^P Ov 

o 

vo 

CN 

CO 

CO 

CN 



ri 

CN 


in 

l O 

C" 

CO 

<Tv 

o 

Cl 

00 

CO 

r- 

i> 

VD 

vo in in 

'N* 

ro 

CO 

rH 

rH 

O 

rH 

rH 

rH 

rH 

o 

O 

O 

o 

o 

o 

o 

O 

rH 

o 

o 

O 

o 

o 

o 

o o o 

o 

o 

o 

O 

O 

O 

O 

o 

o 

o 

oo 

10 

CO 

CN 

in 

0> 

CO 

O 

r- 

in 

00 

m 

00 

00 

CO 

^P CO 00 

CN 

CO 

*p 

OV 

VD 

in 

O 

rH 

CO 

(Ti 

CO 

O'! 

<r 

in 

in 


CO 

O 

in 

<Tv 

rH 

■vP 

o 

ID 

OD rth 

CN 

CO 

(Tv 

in 

CN 

co 

CN 

CO 


O 

rH 

CN 


ID 

C" 

00 

T 

O 

<r 

co 

00 

r 

r- 

ID 

VD 

m in tp 


ro 

CN 

rH 

rH 

o 

rH 

rH 

rH 

H 

o 

O 

o 

o 

o 

o 

o 

H 

o 

o 

O 

o 

o 

O 

o 

o o o 

o 

o 

a 

o 

o 

o 

O 

o 

o 

O 

VO VO 

r 

r- 

co 

r- 

CN 

o 

'P 


cr\ 

CN 

o 

in 

ID 

CO 

O rH rH 


o> 

CN 

o 

rH 

VD 

VD 

in 

o 

VD 


CO 

CO 

co 

ro 

CO 

o 

in 

o 

IP 

CO 

rH 

r~ 

CO 

r- 

in rH in 

<T 

VD 

CO 

in 

CO 

CO 

CN 


CO 

CN 

CN 

CO 

in 

r* 

00 

(T 

o 

<T 

OV 

CO 

t*- 

r- 

ID 

VD 

in 

in in ^p 

CO 

CO 

CN 

rH 

tH 

O 

rH 

rH 

rH 

rH 

O 

o 

o 

o 

o 

O 

rH 

O 

O 

o 

o 

o 

O 

O 

o 

o o o 

o 

O 

O 

o 

o 

O 

o 

o 

O 

o 

LA iH 

CO 

CN 

in 

CO 

o 


in 

oo 

o 

o 

T" 

^P 

O 

V0 

CO ^ CO 

CO 

r- 

VD 

VD 

CO 

o 

CO 

CO 

r- 

T 

rH 

CO 

r- 

rH 

CN 

o 

in 

o 

in 

o 

CO 

in 

CN 

CO 

rH 

<t in o 

VD 

co 

in 

CN 

o 

VD 

o 

rH 

H 


co 

^P 

in 

CO 

cr> 

o 

<T 

<Tv 

co 

CO 

r- 

ID 

ID 

in 

in 

N* 

CO 

co 

CN 

rH 

rH 

O 

rH 

rH 

rH 

rH 

O 

o 

o 

o 

o 

rH 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o o o 

o 

o 

o 

o 

O 

O 

o 

o 

o 

o 

*p Lfl 

r- 


CN 

o 

<T\ 

ro 

rH 

in 

CN 

ro 

o 

r- 

LO 


r- rH CO 

CO 

VD 

CO 

r~ 

m 

ro 

<jp 


VD 

CO 

O 

co 

10 

rH 

o 

CO 

co 


<Tv 


r- 

rH 


CO 

VD 

kP ON t"* 


rH 

CN 

CN 

o 

VD 


o 


rH 

"P 

LD 

r- 

<y\ 

o 

<T 

co 

CO 

r~ 

c~ 

ID 

ID 

in 

in 


^P ^ ro 

CO 

CO 

CN 

rH 

rH 

O 

o 

rH 

rH 

rH 

o 

O 

o 

o 

rH 

O 

o 

o 

o 

o 

O 

O 

o 

o 

o 

o o o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

ro lo 

O 

ro 

o 

r- 

(Tv 

o 

ID 

ID 

in 

CO 

CO 

in 

00 

o 

r~ os vd 

CN 

co 

o 

CO 

in 

CO 

r~ 

o 

rH 

o 

in 

CO 

r- 

o 

CN 


(Tv 


O 

in 

o\ 

CO 

o 

VO 

o 

CO VD rH 

00 

in 

VD 

r- 

in 

rH 

CO 

in 

r- 



vo 

co 

o 

(Ti 

CO 

r~ 

r 

r- 

ID 

in 

in 

in 

'N* 

■p 

CO CO CO 

CN 

CN 

rH 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

rH 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o o o 

o 

O 

O 

o 

o 

O 

o 

o 

o 

o 

cn m 


o 

in 


10 

ro 

^P 

CN 

r- 

CN 

o 

r 

CN 

CO 

oo co r- 

rH 

in 

LO 

rH 


rH 

o 

VD 

VD 


in 

CO 

o 

a\ 

o 

ro 

r" 

CN 

CO 

CN 

ID 

o 

r- 


T" 

C lo rl 

(Tv 

r- 

rH 

CO 

o 

CN 

o 

rH 

CO 

co 

in 

co 

o 

co 

CO 

e'- 

ID 

ID 

in 

in 


^p 

CO 

CO 

CN 

CN CN CN 

rH 

rH 

rH 

O 

o 

O 

o 

O 

o 

o 

o 

o 

rH 

o 

o 

er 

o 

O 

o 

o 

O 

a 

o 

o 

o 

o o o 

o 

o 

O 

O 

o 

O 

1 

o 

O 

o 

o 

cH CO 

o 


<T 

co 

<r 

ID 

ID 

co 

r- 

r 

CO 

ID 

VO 

H 

<T LO O 

o 

o 

VD 

<Tv 


*P 

CO 

rH 

CN 

o 

CN 

o 

VO 


ID 

<T 

co 

CO 

<p 

G\ 

co 

co 

in 

rH 

VD 

id in cn 

o 

a\ 

CO 

CO 


VD 

^p 

CN 

o 

rH 

r- 

o 

00 

r* 

ID 

in 

in 



CO 

CO 

CN 

CN 

CN 

rH 

rl rl rl 

rH 

o 

O 

O 

o 

O 

o 

O 

o 

O 

o 

rH 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o o o 

O 

o 

O 

O 

o 

1 

O 

1 

0 

1 

o 

1 

o 

1 

o 

1 

ro o 

ro 

o 


LO 


ID 

o 

CN 

(Tv 

o 

I — 1 

in 

CN 

rH 

N* OV H 

VD 

in 

CO 

r- 

rH 


(Tl 

ro 

CN 

CN 

o o 


CO 

lO 

m 

CN 

CT> 

ID 

CO 

(Tv 

in 

CN 

Ov 

in 

rH 

<Tv t" rO 

CN 

co 

f" 

o 

CN VD 

VD 

co 

rH 

CN 

O o 

ID 

in 

in 

in 

in 


^P 


CO 

CO 

CO 

CN 

CN 

CN 

rH 1 — 1 rH 

rH 

rH 

O 

o 

O 

o 

o 

o 

O 

O 

O rH 

O 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o o o 

O 

o 

o 

o 

O 

o 

o 

o 

i 

O 

1 

O 

CO 

rH 

CN 

ro 


in 

ID 

C' 

CO 

<Tv 

o 

rH 

CN 

cO 

-Vp 

in vd c-. 

oo 

<Tv 

o 

rH 

CN 

co 


in 

VD 


o 

o 










rH 

rH 

rH 

rH 

rH 

rH rH rH 

rH 

rH 

CN 

CN 

CN 

CN 

CN 

CN 

CN 

CN 


o 


13 



Table 2 Correlation Coefficients Between Like Wind Components (Page 2 of 2) 

INTERLEVEL CORRELATION COEFFICIENTS BETWEEN LIKE WIND COMPONENTS 
ZONAL BELOW DIAGONAL MERIDIONAL ABOVE DIAGONAL 
DATA: KSC SC RAWIN SONDE FEBRUARY 
NUMBER OF PROFILES: 1074 
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Table 3 Correlation Coefficients Between Unlike Wind Components (Page 1 of 2) 
Source: Subroutine QUAD 

Input File: SERCOMKSC.TXT (KSC serially complete Rawinsonde, 1956-74) 
Output File: CORBKSC2.TXT 
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